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ABSTRACT 


It has been shown that good correlation exists between large-scale, high level 
signal returns of bottom acoustic reverberation and bathymetric ridge structures. The 
ultimate goal of this type of analysis would be the removal of propagation effects 
resulting in a large-scale mapping of scattering strengths. Furthermore, analysis with a 
quasi-CW propagation model suggests a strong correlation between the small-scale 
fluctuations in the reverberation signal and the bottom acoustic ensoniflcation, thereby 
suggesting that uncertainties in the predictions of forward propagation may limit the 
resolution of such mapping. By employing broadband modeling techniques, a valid 
representation of the complete time domain forward propagation is provided. Diffuse 
secondary bottom interactions appear to affect the resolution of the primary, direct-path 
interaction at ranges beyond a few water depths. Analysis of data recorded by near¬ 
bottom vertical line arrays (VLA’s) confirms the existence of these secondary, multipath 
interactions in the forward propagation. The exact, two-way travel times fi-om all bottom 
interactions are modeled and the influence of the secondary interactions is quantified. 
Possible ramifications for general sonar system performance are discussed. 
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I. INTRODUCTION 


In recent research conducted using data collected under the Acoustic Reverberation 
Special Research Project (ARSRP) sponsored by the Office of Naval Research, it has 
been shown that good correlation exists between large-scale, high level signal returns of 
bottom acoustic reverberation and bathymetric ridge structures. One may then attempt 
to obtain measures of bottom scattering strengths by removing the influence of 
propagation. Analysis with a quasi-CW propagation model, however, suggests a strong 
correlation between the small-scale fluctuations in the reverberation signal and the bottom 
acoustic ensonification. This suggests that uncertainties in the predictions of forward 
propagation may limit the resolution of such an inversion. The limits of the resolution 
of such a map are also imposed, in part, by the processing and the knowledge of the 
environment. 

By employing broadband modeling techniques, a valid representation of the 
complete time domain forward propagation can be obtained. Despite the use of a 
directional vertical source array, diffuse secondary bottom interactions may affect the 
resolution of the primary, direct-path interaction at ranges beyond a few water depths. 
Analysis of data recorded by near-bottom vertical line arrays (VLA’s) confirms the 
existence of these secondary, multipath interactions in the forward propagation. By 
computing the complex field at the bottom interface from both a vertical line array of 
sources and a point receiver, co-located in a monostatic geometry, the exact, two-way 
travel times and transmission loss from all bottom interactions can be determined. In this 


1 




manner, the influence of the secondary interactions on the measured bottom reverberation 
can be estimated. 

The purpose of this thesis is to examine these effects and attempt to quantify them. 

Other issues are what level of resolution exists within the measured reverberation at 
ranges around one-half of the convergence zone distance and greater, and what level of 
resolution can model predictions hope to provide. Schemes for removing the influence 
of these secondary multipath interactions, if necessary, will also be discussed. Addressing 
these issues will advance the general knowledge of underwater acoustic propagation. 
Quantification of these influences and suggesting schemes to overcome them will lead to 
future improvements in sonar system design. The ability to resolve complex multipaths 
at significant ranges in a covert maimer would provide great advantages to units involved 
in surveillance/power projection operations. 

In Chapter H, some background information on the ARSRP experiment is 
provided. The numerical model used to predict acoustic propagation is also described 
along with other processing techniques used in the data analysis. In Chapter III, the 
existence of the secondary, multiple arrivals is confirmed in the measured data and 
compared with the forward propagation model predictions. The exact, two-way travel 
tim es and transmission loss are then modeled and the influence of the secondary 
interactions is quantified. Findings are summarized in Chapter IV and possible 
ramifications for general sonar system performance are discussed. 
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11. BACKGROUND INFORMATION 


A. THE ARSRP STUDY 

In 1990, the Office of Naval Research (ONR) began an Acoustic Reverberation 
Special Research Program (ARSRP) in the area of underwater acoustic reverberation. The 
scientific plan for the bottom reverberation component of the ARSRP placed special 
emphasis on obtaining a detailed description of the bottom in the "natural laboratory" 
where acoustic experiments were to be conducted [Ref. 1]. The location of the ARSRP 
natural laboratory was selected to be on the west flank of the Mid-Atlantic Ridge 
(approximately 25" to 27" N and 45" to 50" W). Adequate environmental measurements 
could be obtained here in order to understand the connection between sea floor 
characteristics and scattering characteristics. Two geology and geophysics surveys were 
conducted of the region which resulted in data gridded at 200m ^ 200m. More detailed 
analysis was conducted at locations of particular interest providing information on the 
small-scale roughness spectra and bottom composition. 

The ARSRP has conducted two acoustic field experiments in the region. An 
acoustic reconnaissance experiment was carried out in July-August 1991 aboard the 
Research Vessel (R/V) Cory Chouest [Ref. 2]. The major bottom reverberation 
experiment of the ARSRP was carried out in July 1993 [Ref. 3]. Three ships participated 
in the experiment - the R/V Cory Chouest, the R/V Alliance, and the R/V Knorr. The 
Cory deployed a vertical line source array and a horizontal line receive array, the 
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Alliance deployed a flextensional source (vertical pair) and a horizontal line receive array, 
and the Knorr deployed both near-bottom vertical line receive arrays and recording 
packages as well as the; DTAGS (Deeply Towed Acoustics Geophysics System). Fig. 1 
shows the ship tracks for the Cory and Alliance. A major emphasis of the experiment 
was on the collection of bistatic scatter from Sites A, B’, and C’. Both the ship tracks and 
the transmission sequences for the Cory and Alliance were designed so that each ship 
would collect monostatic scatter from their own pings along with bistatic scatter from 
each other. 

During the 1993 acoustics experiment, both the Cory and the Alhance made water 
column temperature and conductivity measurements. The Cory deployed 63 expendable 
bathythermographs (XBTs) as well as four expendable CTDs (1000m). The Alliance 
deployed 39 XBTs and made conventional CTD casts to full ocean depth at five locations. 

B. THE BEAMFORMER 

TTie beamformer developed to process the transmission loss data output from the 
acoustic propagation model was written as a Matlab m-file. The m-file has been written 
to process the data from a 63-element vertical line array, which simulates the near-bottom 
vertical arrays used in the ARSRP study. Inter-element spacing of the simulated array 
is 1.95 meters. The actual arrays used in the study consisted of 59 elements, with an 
inter-element spacing of 2.07 meters. The important consideration for comparison of 
simulated and actual data is that their positions and total lengths are in agreement. The 
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Figure 1. Ship tracks of R/V’s Cory Chouest and Alliance during 1993 main acoustics experiment. From Ref. [4]. 























data is processed to produce a vertical arrival-angle versus time magnitude plot of the 
received sonar signal at the vertical array. 

The beamformer is written for the specific output of the propagation model’s 
broadband simulations. This data corresponds to a frequency bandwidth of 66.66 Hz with 
a center frequency of 227.5 Hz and a total of 1024 frequencies. Thus, incremental 
frequency spacing is 66.66/1024 = 0.065104 Hz. This results in a total Fourier- 
transformed time window of 15.36 seconds and a time resolution of 0.015 seconds. 

These 1024 frequencies are placed in the positive portion of a frequency window 
which extends from -533.33 Hz to +533.33 Hz, with all other frequencies zero padded. 
A Blackman window is applied to the non-zero frequency information, and then the entire 
frequency window is transformed to the time domain resulting in a time window of 15.36 
seconds long, but a time resolution of 9.4 x 10"^ seconds. 

This time window is broken into thirty-one overlapping windows to take advantage 
of the increased time resolution and minimize sidelobe effects m time. Each overlapping 
window consists of 1024 sampled times, which are transformed back to the frequency 
domain where the frequency resolution is now 1.043 Hz. A Blackman window is now 
applied to the 63-element spatial array and transformed to the vertical wavenumber 
domain. Arrival angle information corresponds to wavenumber information as follows 
[Ref. 5]: 


m = 


(N+Z)dcost|f 

m 

X 


(N+Z)dcos\|r f 

m 

C 


( 1 ) 
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where m is the wavenumber bin, N is the number of array elements, Z is the number of 
zeros that the FFT is padded with, d is the inter-element spacing, i|rn, is the arrival angle 
relative to endfire, A is the wavelength, f is the frequency, and c^, is the reference sound 
speed of the water. 

Once the arrival angle information is determined, the frequency window is finally 
transformed back to the time domain. The time domain information in the overlapping 
portions of the time windows are incoherently summed and the final matrix of arrival 
angle versus time is converted to transmission loss in dB (unnormalized) for plotting. 

C. THE ACOUSTIC PROPAGATION MODEL 


The full-wave propagation method used to model the acoustic propagation for 
these studies uses a parabolic equation. Since the introduction of the parabolic equation 
method into the field of underwater acoustics in the 1970’s by Hardin and Tappert [Ref. 
6], the PE method has become a popular wave-theory technique for solving range- 
dependent propagation problems in ocean acoustics. 

1. Standard PE Derivation 

The starting point in the derivation of the standard parabolic equation is the free- 
field Helmholtz equation for a constant-density medium in cylindrical coordinates (r,(|),z) 
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and for a harmonic field of time dependence expf-iwt). 


a^p lap a^p 2 2 „ 

— +-+- +k n p= 0, 

2 r ar 2 0 ^ ’ 


ar 


dz‘ 


(2) 


where azimuthal symmetry is assumed, p(r,z) is the acoustic pressure, = co/Cq is a 
reference wavenumber, n(r,z) = C(/c(r,z) is the index of refraction, and c(r,z) is the sound 
speed describing the environment. 

Assume the solution to Eq.(2) takes the form 


p(r,z)= r), 


(3) 


which is an outgoing cylindrical wave solution. The envelope function i}r(r,z) is assumed 
to be slowly varying in range. The Hankel function, which satisfies the Bessel 
differential equation 


a^iy^”(rr) . dH^'\k r) 

0 0^ 1 0^0'^ 


(1). 


dr^ 


dr 


^klH^\kx)= 0 , 

0 0 0 


(4) 


can be replaced by its asymptotic form for )fe(,r » 1, 



% r) = 




Tzk r 
0 


i(*^r-ir/4) 

e 


(5) 
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Substituting Eq.(3) into Eq.(2), and using the Hankel-function property of Eq.(4), gives 




H^^\k r) 
0 ^ 0 


ar 


1^ 
+ — 

ry 


aiif a^ijt , 2 2 

— + —-+k (n - 1 ) 1 ]; = 0. (6) 

^ dz^ ® 


Making the far-field assumption, » 1, Eq.(5) can be used to obtain the simplified 
elliptic wave equation 


a^ijt 

ar^ 


atjt 

+ Tik — + 

0 ar 


a^i|t 

dz^ 




(7) 


Next, the paraxial approximation is used to arrive at the standard parabolic wave 
equation. This small-angle approximation is expressed by 


a^t|j 

ar^ 


aiit 

clik —. 

0 ar 


( 8 ) 


A qualitative justification for this approximation can be made by noting that the main 
radial dependence of the field is contained in the Hankel function through the term 
exp(i^), while the envelope tjr will vary slowly with range over a wavelength A. This 
can be expressed mathematically as 0t|f/ar < ~ ijr, or a^x|f/ar^< ik^dT}f/dr , and 

consequently the paraxial approximation is justified. 
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Making use of Eq.(8) in Eq.(7) results in 


2ik 


0ilr 

di 


dz^ 


* (n"-l)l|l = 0, 


or 


(9) 


3t|r 


i d^tjj 


ik 


which is the standard parabolic equation. 


2. The Split-Step Fourier Algorithm 


A second-order differential equation in r has now been reduced to a first-order one 
in r, allowing solutions by way of a non-iterative marching algorithm. Eq.(9) may be re¬ 
written as 


i 0ijr 

k di 
0 



t = (r 


op 


op 


( 10 ) 


where the operators 


T 

op 



i d 


\2 


k dzj 
0 


1 

2k^ dz^ 
0 


( 11 ) 
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and 


^op " = --i (n^ -1). (12) 

This representation of the operators as kinetic and potential energy operators is 
especially insightful when it is desired to form the ray equations which have Hamiltonian 
form [Ref. 7]. In Eq.(lO), the function t|; is a vector (in z) in Hilbert space. The values 
of tj; can now be expressed as 

i|f(r + Ar) = $(r)tlr(r). (13) 

For all propagation predictions presented, we have employed the University of 
Miami Parabolic Equation (UMPE) model of Smith and Tappert [Ref. 8]. The UMPE 
model uses the spUt-step Fourier (SSF) method [Ref. 6] to compute the PE solution. This 
is accomplished by approximating the propagator function ^(r) by 


where 


$(r) » e 


-ik H (r)Ar 

0 op^ ' 


H 

op 



1 

Ar 


I 


r + Ar 




op 


(r') 


(14) 


(15) 


This can be approximated by 


H 

op 



(r + ^Ar), 


(16) 


11 




or simply 


(17) 

which are referred to as "centered" and "endpoint" schemes, respectively. Presumably if 
the range step Ar is small enough the differences between the solutions are negligible. 
The UMPE model uses the approximation of Eq.(17). 

The operator U^p is a multiplication operator in z-space and is therefore a diagonal 
matrix. The operator T^p is not diagonal in z-space, but is diagonal in wavenumber space. 
By separation of H^p into the two components T^p and U^p the propagator function 
becomes 


0(r) = e 


-ikAtT -ik AiU 
0 ^ ^ 


where the Baker-Campbell-Hausdorff expansion [Ref. 9] 


(18) 


e ^ = e ^ ^ ^ + • • • (19) 

has been used. Since both T^p and U^p are small, their products are of second order and 
assumed negligible. 

The general algorithm behind the PE/SSF implementation is then as follows. The 
PE field function i|f is specified at some range r in the z-domain. A transformation is 
made to the A:-domain followed by a multiplication of the k-space operator e . 

The result is then transformed again to the z-domain and is followed by a multiplication 
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-ik LtU 

of the z-space operator e . The final result is the field function at r + Ar. The 

FFT subroutine employed in the numerical code assumes the convention 


i|r(z) = FFT(tIr(^)) = 


( 20 ) 


and 


m = 


PFT(i|f *(z)) * = I* dz. 


Therefore, the PE/SSF implementation can be represented by 


( 21 ) 


\|f (r + Ar) = e 


-ik AtU (t,z) 
0 


^ FFlje 




X [FFr(i|(‘(r,z))]‘ 


( 22 ) 


where, in fc-space, the operator in Eq. (11) is given by 


t (k)=ll^]\ 

op l^k' 


(23) 


and Uop is as defined in Eq. (12). 

The UMPE model allows the user to choose from five different sets of operator 
forms. For the modeling used in this thesis, the higher order wide-angle approximation 
(WAPE) form [Ref. 9] has been used which is accurate to greater than ± 20° from the 
horizontal. The operator forms for the wide-angle approximation are 

dz^ ^ 


WAPE 


0 

2 
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and 


U 

WAPE 


= -(n-1). 


In wavenumber space Eq.(24) may be written as 


(25) 




-1 


= 1 - 


Note that modes with k > are evanescent since 





(26) 


(27) 


In the modeling done for this thesis, measured sound speed profiles from the 
Alliance CTD casts are used as input. The bottom is typically modeled as a single 
interface at depths extracted from the gridded database. The acoustic bottom properties 
are chosen to be representative of this area and have no lateral variability due to lack of 
geophysical data. Various analytical formulations of source functions are available; the 
point source and vertical line array were used in this work. Broadband UMPE model 
runs were performed to allow for more accurate investigations of the time domain 
response of the medium. 
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3. Computing One-Way (Forward) Travel Times 


The UMPE code calculates the time domain arrival structure at a given range in 
a straight forward manner. The time harmonic acoustic field assumed by the model is 
given by 


P(r,z,(»)t) =p (r,z)e 

(0 


(28) 


The representation in the time domain is then simply 


P(r,z,t) =FFt[p (r,z)l = fp (r,z)e"'“^do), 

, u) J J <0 


(29) 


which follows from the realization that Eq.(28) is a single firequency component of the 
full-spectrum time dependent field given by Eq.(29). 

To compute the arrival structure at some fixed range r=R, the model computes 
and stores the complex field values Pujr(z) for many frequencies and then Fourier 
synthesizes them to obtain P (z,t), the set of complex pressure values in time/depth 
space. Because the EFT assumes inputs over the frequency band 4 - BW/2 to f^ + BW/2, 
where BW is the bandwidth of the acoustic source and f^ is the center frequency, high 
frequency calculations can become computationally burdensome even for small 
bandwidth fields. Therefore, the UMPE model frequency shifts the center frequency to 
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d.c. The calculation becomes 


_ / r -j(u-w )t 

Pjfet) = I p^^(z) e • d(0, (30) 

where o)o = 27tfo. 

In practice, the model computes the FFT of (1 /^) ^^(z). Because this neglects 

ik R ) 

the overall phase factor g ° = e *, the time domain is heterodyned around the value 
to = R/Co- Arrival times are then given as values of "reduced time", or (t - O. 

To minimize the effects of sidelobes due to a flat frequency spectrum, and of 
aliasing or wrap-around due to the finite BW, the UMPE model creates a Hamming 
window which is applied to the entire bandwidth. This method of tapering greatly 
reduces the effect of side lobes. Wrap-around effects may still be evident if the time 
window is not wide enough. The total length of the time window is 


T = 


N 

BW 


(31) 


where BW is the bandwidth of the source and N is the number of frequencies for which 
the model computes solutions, or the frequency FFT size if the FFT is padded with zeros. 
The frequency resolution is 


Ao) 


BW _ 1 

~W~T 


(32) 


so the wider the time window needed, the finer the bandwidth must be resolved and the 
more total runs at different frequencies must be performed. 
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4. Computing Exact, Two-Way Travel Times 

The time domain arrival structure of the acoustic field at a given point in range 
and depth can be constructed from the Fourier synthesis of the full-spectrum complex 
field values calculated by UMPE at that range and depth over the frequency bandwidth. 

In this manner, for every point on the bottom, a time series can be constructed 
which would show the arrival structure of the acoustic field from some source at a 
distance r = R from the point on the bottom. 

Imagine a second source co-located at range R from the point on the bottom. If 
we label the bottom point b, and the two co-located sources r and s, then the arrangement 
would be as shown in Fig. 2. 


S 



Figure 2. Source and Receiver Arrangement 
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The time series at a bottom point b from a source at location s could be 
constructed as previously shown. Also, the time series at receiver location r from a 
source at location b could be constructed in the same manner. By invoking reciprocity, 
the complex field at location r from a point source at b over a given frequency spectrum 
is the same as the complex field at location b from the same point source at r. 

ITie time domain arrival stracture of the complex acoustic field at location r from 
a source at point s would be a convolution of the time-domain arrival structures from s 
to b and from b to r: 


= /ilr^^(t-T,R) ijr Jt,R) dx (33) 


or, in the frequency domain: 


t^Cco.R) = iIr^/w,R) X tjr Ja),R) 


(34) 


where the subscript "2" stands for two-way, and ((«),R) = tjr ^(a),R)by reciprocity, 
tir (&),R)and tk (<«),R) are the complex acoustic fields output by the UMPE model for 
all frequencies in the spectrum. 

This development has assumed reflection of tbe field at only one bottom point. 
In reality, the arrival structure at point r is the summation of the convolved arrival 
structures for all range steps: 


2Ar , 4Ar 

-,r+Ar) + il; 2 (t+-,r+2Ar) + ... 

c c 

0 o 


( 35 ) 
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where all time series developed for each step in range must be summed. For any given 
range R, the center time from the Fourier synthesis will be T,. = R/Co with a time window 
width of AT = l/6f where 6f is the incremental frequency spacing between each CW 
calculation. A range step was used to produce changes in equal to an integer number 
of 6t = 1/BW, the corresponding incremental sampling of the time window, i.e. 



In order to develop the full time series, the exponential factor immediately above 
will be left out of the actual calculation, and each subsequent range step contribution will 
simply be shifted in time so that its center time is t = 2Ar/c. 
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III. DATA ANALYSIS AND MODEL PREDICTIONS 


A. THE MODELED STUDIES CONDUCTED 


In order to perform model predictions which can be compared to actual five 
segments were chosen from the ARSRP study [Ref. 3] which correspond to five distinct 
arrangements of towed source transmitting to a near-bottom vertical array. For all cases, 
the model source was placed at 181 meters depth, a source level of 232 dB re pPa/m was 
used, and bathymetry was extracted along the bearing from source to receiver in order to 
accurately simulate the true environment. The sound speed profile in the water column 
used by the model was taken from one of the CTD casts in the region, and the bottom 
was assumed to be homogenous with a sound speed, density, and attenuation of 1800 m/s, 
2.0 g/cm^, and 0.2 dB/m/kHz respectively. 

For Segments 016,031, and 046, the near-bottom vertical line array (VLA) was 
stationary as the Cory’s source moved around it. The UTM coordinates (zone 23, based 
on the WGS-84 ellipsoid) for the near-bottom array were (377.43, 2888.51). The source 
for Segment 016 was located at UTM coordinates (383.82, 2901.80) so that the range 
from source to receiver was 14.75 km. The source for Segment 031 was located at 
(398.03, 2899.89) and the separation range was 23.53 km. The source for Segment 046 
was located at (388.66, 2885.01) and the separation range was 11.76 km. 

For Segments 544 and 564, the near-bottom VLA was located at UTM coordinates 
(192.27, 2944.86). The source for Segment 544 was located at (223.2, 2953.54) and the 
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separation range was 32.12 km. The source for Segment 564 was located at (210.36, 
2932,48) and the separation range was 21.92 km. 

Plots of the transmission loss fields of the ocean cross-sections along the radials 
of the five segments were made at the center frequency of 227.5 Hz. These fields are 
calculated by the UMPE model and show the propagation of the acoustic energy from the 
source outward. 

Segment 016’s arrangement, shown in Fig. 3, places the near-bottom array at 14.75 
km, which is prior to one-half convergence zone (CZ) distance. The bottom profile for 
this radial is relatively smooth with a slight upslope along the latter half of the 
calculation distance. 

Segment 03rs arrangement, shown in Fig. 4, places the near-bottom array at 23.53 
km, just prior to one-half CZ distance. The propagation in this case is bottom-limited at 
25 km. The bottom profile for this radial shows an almost flat bottom along the mid¬ 
point with gradual rises of approximately 1500 meters near source and receiver locations. 
The forward scatter resulting from interaction with the bottom can be dimly seen in the 
top right of the figure of which some reflects off the surface back down to the location 
of the array near the bottom. 

Segment 046’s arrangement, shown in Fig. 5, places the near-bottom array at 
11.76 km which is well before the one-half CZ range. The bottom profile shows several 
larger-scale bottom features in the first half of the radial distance. 

Segment 564’s arrangement, shown in Fig. 6, places the receiver array at 21.92 
km, prior to one-half CZ distance. The bottom profile shows a gradual downslope from 
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Figure 3. Segment 016 Transmission Loss Field for the Array Source. 
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Figure 4. Segment 031 Transmission Loss Field for the Array Somce. 
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Figure 5. Segment 046 Transmission Loss Field for the Array Source. 
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Figure 6. Segment 564 Transmission Loss Field for the Array Source. 
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5-10 km and a gradual upslope over the last quarter of the radial distance, with no large- 
scale features interrupting the relatively smooth profile. 

Segment 544’s arrangement, shown in Fig. 7, has one of the most interesting 
bottom profiles with several major large-scale features from which reflected energy paths 
can be clearly seen. The receiver array is located at 32.12 km, just past one-half CZ 
distance. In this case propagation is bottom-limited at -30.5 km, prior to the array 
location. Figure 8 shows the UMPE model prediction of the approximately 200-255 Hz 
broadband pulse arrival transmission loss structure over the entire water column at the 
range of the near-bottom VLA. The separation between primary and secondary bottom 
interactions is observed to be roughly 1.5 seconds. 

B. EXISTENCE OF SECONDARY MULTIPLE ARRIVALS 

Imtially, the transmitted pulse in all cases can be considered a single wavefront 
in time. The highest energy level is transmitted horizontally from the source and occupies 
the higher, more shallow end of the pulse. As the pulse propagates the bottom end 
undergoes a bottom reflection which begins to spread the signal’s range resolution. 
Further dowmrange, this forward scattered energy reflects from the surface and again 
interacts with the bottom. This propagation is explained in detail in Ref. [4]. 

The range resolution of this secondary interaction has been much degraded as the 
pulse has become spread out over some range. Also, this secondary interaction is 
occurring at the same time the primary direct-path pulse is ensonifying a patch of seafloor 
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further downrange. Since the reciprocal paths of both the primary and secondary 
interactions back to a monostatic receiver platform will arrive at the same time, discrete 
scattering events at the position of a primary interaction will be contaminated or 
combined with scattering events at different places on the seafloor. These secondary 
interactions occur at higher angles of incidence and, because of this, may have their 
scattering strengths enhanced. 

The effects of these secondary interactions are most significant at ranges of one- 
half CZ and greater. In cases where the propagation is not bottom-limited, or where the 
primary portion of the pulse is forward scattered, the highest level of ensonification 
beyond the range of primary interactions will be due to the secondary interactions instead 
of the primary. The returns from these interactions can result in reverberation features 
being mapped beyond the true point of contact by several kilometers. 

By looking at the data recorded by the near-bottom VLA’s in the vicinity of one- 
half CZ, such phenomena can be verified. The ability to beamfotm the data allows the 
primary, near-grazing interaction to be distinguished from the secondary, higher-angle 
arrivals undergoing multiple reflections. 

Data from the five simulated segment runs has been extracted over the aperture 
of the actual near-bottom VLA (120 meters long, approximately 100 meters off of the 
bottom). Vertical spacing of the model data required 63 depth elements to be extracted 
to simulate the complete aperture of the near-bottom VLA. This data was then 
beamformed as explained in Chapter II to yield plots of transmission loss arrival 
structures for the five segments. 
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Figure 9 shows the beamformed predicted arrival structure for Segment 016. 
Values of transmission loss in dB (unnormalized) outside the plotted scale are assigned 
the maximiun or minimum values for the plot as appropriate. The exact normalization 
factors for transmission loss (TL) values are unknown due to the process of beamforming 
the data, but all modeled and beamformed values are accurate relative to each other. 
Positive angles indicate propagation from the bottom towards the surface, i.e. energy 
coming up from below these element depths. The initial direct path arrival is so intense 
that the sidelobe structure has spread over multiple adjacent angles and times. The next 
set of arrivals occurs roughly 2.5 seconds later and exhibits a symmetry about grazing 
with lower intensity and more diffuse arrivals from below (positive angles). The arrivals 
from below are more diffuse due to an extra reflection from the bottom just prior to 
reaching the array. These secondary arrivals tend toward higher arrival angles with 
increasing time as expected due to the increased path lengths they must travel. In 
addition to these features, there is also some structure which tends toward more grazing 
angles with time. This seems to originate from a point of high-intensity incidence which 
also produces the expected increasing angle with time features. The cause of this 
unexpected feature is not well understood, though all features must be the result of 
forward scatter since the model does not compute any back-scattered effects. 

The match-filtered and time-domain beamformed LFM signal arrival structure 
received by the VLA for Segment 016 is shown in Fig. 10. The primary arrival occurs 
at 0.4 seconds and shows earliest arrivals from above (negative angles) with later arrivals 
occurring as the arrival angle passes through zero to positive values, which also matches 
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with model predictions. The structure seen at 10* and 1.3 seconds is most likely 
backscatter from some bottom feature just past the VLA. This structure does not show 
up in Fig. 9 as the model only calculates forward scatter. Secondary arrivals occur 
approximately 2.8 seconds after the primary and exhibit the same structure as described 
above for the model predictions of Fig. 9. 

Figure 11 shows the beamformed predicted arrival structure for Segment 031. The 
general pattern is the same, though the second set of arrivals, which occurs roughly 1.5 
seconds after the primary arrival, is more spread out in time as would be expected for the 
longer propagation range. 

The match-filtered and time-domain beamformed LFM signal arrival structure 
received by the VLA for Segment 031 is shown in Fig. 12. The primary arrival occurs 
at 0.25 seconds, there exists some backscatter structure at 1.5 seconds, and secondary 
arrivals show up 1.8 seconds after the primary arrivals. This shows good overall 
agreement with the model predictions shown in Fig. 11. 

Figure 13 shows the beamformed predicted arrival structure for Segment 046. 
Secondary arrivals first appear 2.8 seconds after the primary. The match-filtered and 
time-domain beamformed LFM signal arrival structure received by the VLA for Segment 
046 is shown in Fig. 14. Again, arrival angles and time between primary and secondary 
arrivals agree very well. 

Figures 15 and 16 show the beamformed predicted arrival structme for Segment 
544. No actual beamformed data was available for this segment. The model predictions 
conducted have a frequency spacing which results in a total time window of 15.36 
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Figure 11. Segment 031 Beamformed Arrival Stmcture. 
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Figure 12. Segment 031 Beamformed VLA Arrival Stmcture. 
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Figure 15. Segment 544 Beamformed Arrival Structure. 
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Figure 16. Segment 544 Complete Beamformed Arrival Structure. 
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seconds, and this longer time window is shown in Fig. 16. There exists a third, fourth, 
and possibly a fifth distinct arrival structure, each of which arrives at progressively higher 
angles and with less time separating each subsequent set of arrivals. Comparison of the 
full time windows of beamfotmed modeled data shows that in general, the farther out in 
range from the source, the shorter the time interval between the primary and first 
secondary arrivals; the more spread out in time these secondary arrivals; and the shorter 
the time between subsequent sets of secondary arrivals, all of which must result from 
forward scatter. For the purpose of this thesis, only the secondary multi-path arrivals 
which occur at or close to the same time as the primary arrival are of concern. 

Figure 17 shows the beamformed predicted arrival structure for Segment 564 and 
Fig. 18 shows the match-filtered and beamformed data from the near-bottom VLA for this 
segment. Given that the predicted arrival structures for all five segments are the result 
of only forward propagation calculations, there is sigmficant agreement between predicted 
and actual arrival structures for all cases, showing that these secondary arrivals do indeed 
exist and can affect two-way mapping resolution. 

C. INFLUENCE OF SECONDARY INTERACTIONS 

The existence of secondary arrivals has been clearly shown in the previous section. 
To what extent do these secondary arrivals degrade the return signal at the receiver in a 
monostatic arrangement? 
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In an attempt to answer this, exact two-way travel time arrival signal levels were 
calculated for Segments 016 and 544 as explained in Chapter 11. In order to ensure that 
a Fourier-transformed time window would be large enough to capture the complete 
structure of the two-way propagation out to the farthest range needed, frequency spacing 
for the UMPE modeling runs to determine exact two-way travel times was halved so that 
there were 2048 frequencies centered at 227.5 Hz with a total bandwidth of 66.66 Hz, 
resulting in a total Fourier-transformed time window of 30.72 seconds. 

The exact two-way travel time structures and corresponding bottom profiles for 
Segments 016 and 544 are shown in Figs. 19-22. There is much more structure in both 
Figs. 19 and 21 than exists for the corresponding bottom profiles. TTus is an indication 
of the high degree of complexity in the ensonifymg field. The fact that such complexity 
can exist distinct from any fine-scale bathymetric roughness is described in more detail 
in Ref. [4]. Large features can be discerned by large-scale arrival structures and trends 
in the return signal, but small-scale features correspondmg to the "spilg^-ness" of the 
return signal simply do not exist in the modeled bathymetry. Lobe structure from the 
source array as seen m Fig. 3 is the cause of the return signal’s general appearance of 
very low return up to approximately 4 km, stronger returns from 4 to 10 km, and 
strongest returns from 10 to 15 km, the maximum range considered. 

In an attempt to remove the majority of secondary interactions from the 
propagation, the UMPE code was modified for each segment to effectively remove the 
bottom interface up to a range of five kilometers in front of the near-bottom VLA’s. This 
was done by matching water and bottom densities and sound speeds along the radial 
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range. At five kilometers prior to the VLA, the bottom properties were linearly changed 
so that within three kilometers of the VLA, the bottom properties were the normal values 
used in the previous runs as stated earlier in the chapter. Though the bottom interface 
was removed, energy entering the ocean bottom was still attenuated as before. The 
results of altering the bottom properties for these segments on the propagation of acoustic 
energy can be seen in Figs. 23 and 24 which show the transmission loss fields of the 
ocean cross-sections along the radials of the two segments at the center frequency of 
227.5 Hz. 

Figure 23 shows the propagation along Segment 016. Surface reflections still 
occur, but any energy reflected from the surface near the source would hit the bottom 
interface prior to the VLA in the altered bottom region and be attenuated. With multipath 
energy removed, the propagation paths of the primary energy close to the surface can now 
be clearly seen. The three primary propagating lobes of the vertical array source are 
distinct, and the interference pattern resulting from Lloyd’s mirror effect of the surface- 
reflected energy within the lobes is also more distinct. Figure 24 shows the propagation 
along Segment 544. Bottom-reflected energy just prior to the VLA can be seen to reflect 
from the surface and travel back down to ensonify the region beyond 31 km range. This 
region sees little of the primary energy which is bottom-limited around 31 km. 

Figures 25 and 26 compare the unaltered bottom signal return and the altered 
bottom signal return for Segment 016. At ranges approaching 15 km, the same large 
structure exists in the returns, while in Fig. 25 it can be seen that multipath effects tend 
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Figure 23. Segment 016 Transmission Loss Profile - Altered Bottom. 
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Figure 24. Segment 544 Transmission Loss Profile - Altered Bottom. 
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to fill in the areas of low primary signal return, adding up to 5 dB of return in some 
cases. While the effect of the secondary, multipath arrivals is not very significant at these 
short ranges, it is important to note the considerable differences in details in the region 
of the altered bottom. The distinction can be considered the difference between 
computing the total (incident, reflected, transmitted, and forward scattered) field at the 
interface versus computing only the incident field. To correctly predict reverberation 
levels, the total field at the interface must be computed. 

Figures 27 and 28 compare the unaltered bottom signal return and the altered 
bottom signal return for Segment 544. Again, the general effects of source array lobe 
structure can be seen when compared to Fig. 7. In the range of approximately 17 to 21 
km, signal returns in the unaltered bottom case are as much as 25 dB higher than the 
returns caused by primary interactions alone. By referring back to the CW field in Fig. 
7, it is noted that a ridge structure is near 16 km. This would produce a shadow region 
in the 17-20 km range if only direct path ensonification were considered. However, 
secondary multipaths are predicted to significantly increase the levels in this region. 
Another interesting feature is that the primary return in the unaltered bottom case at the 
range of approximately 21 to 23 km is actually lessened by the multipath effects. This 
is presumably due to the difference between computing the total field versus only the 
incident field. Once again, the lower primary return at the range of 25 km due to a 
shadow zone behind a ridge at approximately 23 km has had as much as 30 dB added to 
it in Fig. 27 from the returns of secondary, multipath interactions being received at the 
same time as the primary return. 
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Figure 28. Segment 544 Two-Way Exact Travel Time Structure - Altered Bottom. 




































Figures 29 and 30 show a magnified view of the imaltered bottom signal return 
and the altered bottom signal return, respectively, for Segment 544 for the ranges of 25 
to 33 km. The secondary signal returns in Fig. 29 could result in some features being 
mapped to incorrect seafloor locations. The structure in Fig. 29 from 31 km to 32 km 
which exhibits the gradual upslope then downslope is actually from secondary returns 
occurring approximately 1.5 seconds later (see Figs. 8 and 15) from the large-scale ridge 
face at 30 km and would result in mapping another ridge structure at this location. This 
secondary "echo" was previously predicted in Ref. [4]. 
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IV. SUMMARY 


It has been shown that the secondary arrivals predicted by the UMPE model do 
exist in the measured data. The exact two-way travel time calculations show that the 
levels of ensonification resulting from these secondary multipaths can be significant 
compared to the primary direct-path interaction at ranges near one-half GZ and greater. 
Furthermore, the exact two-way levels display a highly complex nature which are due to 
complicated ensonification structures rather than small-scale bathymetric features. The 
overall implication of this is that there are limitations on the ability to resolve discrete 
structures at ranges beyond a few water depths. These limitations arise from system 
limitations, processing techniques, and the inherently complex nature of the propagation. 

One way to overcome some of these limitations in mapping resolution might be 
to use a vertical receive array, or some form of volumetric array, so that the ability to 
beamform in the vertical exists. Most secondary multipath interactions, at least those 
occurring during the return of the acoustic signal, will be arriving at higher angles than 
the direct-path returns so that they can be selectively rejected from the signal return. 

Another aspect of the propagation which may affect the resolution of such a 
system is coupling between radials. This has not been considered in this research. 
Further studies incorporating azimuthal coupling should be conducted in an attempt to 
better understand the inherent problems associated with attempting to precisely map 
bottom bathymetry in this manner. 
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Finally, the shallow water environment is much different from the deep ocean 
environment and the modeling done for this thesis could be applied to shallow water 
cases to determine the effects of these secondary interactions. In a shallow environment, 
most of the propagation will be bottom-limited so multipaths will be the primary source 
of bottom ensonification. With the shift in emphasis toward the littoral environment for 
the Navy of the future, these issues will need to be addressed for unproved sonar system 
performance on the submarines and ships of tomorrow’s Navy. 
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